Spontaneous parity breaking of graphene in the quantum Hall regime 
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We propose that the inversion symmetry of the graphene honeycomb lattice is spontaneously 
broken via a magnetic field dependent Peierls distortion. This leads to valley splitting of the n — 
Landau level but not of the other Landau levels. Compared to quantum Hall valley ferromagnetism 
recently discussed in the literature, lattice distortion provides an alternative explanation to all the 
currently observed quantum Hall plateaus in graphene. 



Recent experiments have revealed peculiar quantum 
Hall (QH) effects in graphene, a single atomic layer of 
graphite [J, Q. These measurements are understood as 
single electron effects and unusual QH features can be 
traced back to the relativistic-like dispersion relation of 
electrons in graphene and to their twofold valley degen- 
eracy. In particular, the observed plateaus in the Hall 
conductivity at filling factor v = ±2:±6; ±10 can be eas- 
ily understood in this framework [3]. Following these 
pioneering experiments, Zhang et at 0] discovered 
new QH plateaus at v — 0;±1;±4, which several au- 
thors 0, 0, 0, H @ attribute to valley (and spin) fer- 
romagnetism, relying on interactions between electrons. 
However, the absence of plateaus at v = ±3; ±5 is in- 
triguing in this respect and cast doubts on this interpre- 
tation. Alicea and Fisher [l| propose that disorder might 
be so strong in current graphene samples as to destroy 
exchange interactions and therefore ferromagnetism [a]. 
If this is indeed the case, one still has to explain the ori- 
gin of the extra plateaus. Alicea and Fisher suggest a 
valley splitting mechanism relying essentially on lattice 
scale electron repulsion while neglecting exchange inter- 
actions in a dirty graphene sample [9|. Another scenario 
relying on long-range electron interactions and leading to 
an excitonic valley gap is the so-called "magnetic cataly- 
sis" [13] ■ In the present paper, we take a different route 
and assume from the outset that direct interactions be- 
tween electrons do not play a major role. We explore the 
possibility that all the plateaus observed so far could be 
understood as integer QH states resulting from electron- 
lattice interaction effects. The new input of our model 
is a magnetic field driven out-of-plane lattice distortion 
lifting the valley degeneracy. 

Graphene is a honeycomb lattice of carbon atoms: a 
two dimensional triangular Bravais lattice with a basis of 
two atoms, usually referred to as A and B. The distance 
between nearest neighbor atoms is a = 0.14 nm and the 
lattice constant is a\fi. Experimentally, graphene sheets 
of area A ~ (3 — 10 /zm) 2 are deposited on Si02/Si sub- 
strate. Applying a gate potential V g via the substrate 
allows one to control the electronic filling of the graphene 
bands. The number of induced electronic charges is given 
by N c — VgCg/e where the capacitance per unit area can 
be estimated as C g /A ~ e r eo/d 1.2x 10 -4 F/m 2 , where 



— e < is the electron charge, e r w 4 is the silicon oxide 
dielectric constant and the thickness d ~ 300 nm [l|, l2j. 

In order to study the electronic properties of graphene, 
we use a standard nearest neighbor tight-binding model 
11 1 with hopping amplitude t w 3 eV It describes 
the hopping of electrons between 2p z carbon orbitals. 
There is one electron per carbon atom. If we call N p the 
number of plaquettes (or unit cells) , there are 2N p elec- 
trons in the sample under zero gate potential. The first 
Brillouin zone is hexagonal and of its six corners, only 
two are inequivalent and usually called K and K' . We 
choose K = 47r/(3v / 3a)u 2; and K' = — K. The resulting 
band structure features the merging of the conduction 
and valence band at precisely these two points: graphene 
is a two valley (K and K') zero-gap semiconductor. 
Near these so-called Dirac points, the electrons behave 
as charged massless Weyl (or chiral Dirac) fermions with 
Fermi velocity vf — 3at/2h w 10 6 m/s playing the role 
of an effective light velocity in the relativistic-like disper- 
sion relation Sk = ±top|k|. When the gate voltage is 
zero, the "big band" (valence plus conduction band) is 
half-filled: the Fermi level is right at the Dirac points. 

Adding a weak perpendicular magnetic field B± , such 
that the flux per plaquette is much smaller than the 
flux quantum (f>o = h/e, McClure first obtained the rel- 
ativists Landau levels (LL) of graphene see also 
Ref. [3, flj| . Including the Zeeman effect, the LL in the 
Dirac equation approximation read 



sgn(n) 
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where the "cyclotron energy" is hiu c = v2Hvf/Ib, the 
LL index n is an integer, the spin projection along the 
magnetic field axis is a = ±1, the Bohr magneton is 
jiB = eh/2m, with m the bare electron mass, and the 
effective g-factor is g* ~ 2 close to its bare value Q • The 
magnetic length is defined as usual by Ib — \/h/eB±. 
If the Zeeman splitting is negligible, each LL has de- 
generacy 4iV0. The total number of flux quanta across 
the sample — B±A/<fio gives the orbital degeneracy. 
The factor 4 accounts for spin 1/2 and twofold valley de- 
generacy. We call v — Nc/N^ — CgVg/eN^ the filling 
factor. When the gate voltage is zero, v = and the 
n = (central) Landau level (CLL) is half-filled as a re- 
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suit of particle-hole symmetry leading to the remarkable 
fact that the number of electrons in the CLL is 2N^. 

We now consider a spontaneous out-of-plane lattice 
distortion that - in presence of a substrate - breaks the 
inversion symmetry of the honeycomb lattice and pro- 
vides a mechanism for lifting the valley degeneracy. As- 
sume that the A (resp. B) sublattice moves away (resp. 
towards) the substrate by a distance 77 [22]. Electrons are 
still described by a honeycomb nearest neighbor tight- 
binding model, however the two atoms in the basis now 
have different on-site energies The energy on atom 
A/B is called ±M following Haldane jl4J, who calculated 
the LL of such a system. Close to the Dirac points, it 
reads 
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where a = ±1 is the valley index corresponding to the 
Dirac points aK. In terms of the low-energy effective the- 
ory, the distortion means that the Weyl fermions sponta- 
neously acquire a finite mass. Note that the on-site en- 
ergy difference lifts valley degeneracy for the CLL only. 
In addition the effect of a nonzero on-site energy M on 
each n ^ LL is very weak, of order M 2 /Hv F eB± ~ 
5.10 -4 for a typical magnetic field ~ 10 T as we will see. 
We could therefore set M = in the n^O LL and use 
the approximate Eq. fTJl instead of Eq. However, 
we shall see below that in order to compute the lattice 
distortion it is important to keep Eq. (J2j) - 

Such a lattice distortion spontaneously occurs because 
it lowers the total energy, in a way similar to Peierls's 
mechanism except for the magnetic field playing an 
essential role here and for the crystal being two rather 
than one dimensional. Assume that the last partially 
filled LL is n = (i.e. the gate voltage V g is such that 
\v\ < 2). We show that in this case it is always favor- 
able to slightly distort the lattice provided there is a per- 
pendicular magnetic field [24j]. The distortion lowers the 



electronic energy. This energy lowering comes both from 
the CLL, which gives an essential contribution, and also 
from all the n < LL, which contribute in a less im- 
portant way as we explain below. There are (2 + v)N,j ) 
electrons in the CLL. They contribute an energy gain 



E n = = -A^(2 - \v\)M 



(4) 



because when v < 0, all (2 + v)N,p electrons gain each 
an energy M and when v > 0, 2N ( j > electrons gain each 
an energy M but the remaining electrons loose each 
an energy M. This energy gain depends on the mag- 
netic field through N^. In addition, the energy gain is 
linear in the out-of-plane distortion 77 because the on-site 
energy is proportional to the distortion, as we discuss 



below: M = D-q, where D is a proportionality constant, 
akin to a deformation potential. The other 2(N p — N^) 
electrons that fill the n < LLs, also contribute to the 
energy lowering. Each of them gains a small energy com- 
pared to what an n — electron gains, as discussed in 
the preceding paragraph, but as there are many more of 
them, about 2(N p — N^) w 2N p , we can not neglect their 



contribution, 
find 



In the Dirac equation approximation, we 
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where the numerical factor 7 = 3 1 / 4 / ~ 0.74 [251 ] . 
This energy gain is quadratic in the distortion, and there- 
fore smaller than E n= o at small distortion, and indepen- 
dent of the magnetic field. Actually, this term represents 
the full electronic energy gain for a lattice distortion un- 
der zero magnetic field. In the end, adding E n< Q to E n —o, 
we see that the larger the magnetic field, the larger the 
electronic energy gain. 

The distortion costs an elastic energy 



-^elastic — NpGf] 
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where the out-of-plane distortion is assumed to be small 
77 <C a and G is an elastic constant. As E n<0 and Elastic 
are both quadratic in the lattice distortion, we introduce 
a renormalized elastic constant G' — G — ^aD 2 /Hvf and 
write an effective elastic energy: 



-^elastic + E n< Q — N p G'r] 2 



(7) 



The effect of the n < electrons is to reduce the lat- 
tice stiffness and therefore to enhance the distortion. We 
take it as an experimental fact that there is no sponta- 
neous out-of-plane distortion in absence of perpendicular 
magnetic field, see also which means that G' > 

We now estimate the two constants D and G. From 
the frequency ojq/2ttc ~ 800 cm -1 of the graphite out-of- 
plane optical phonon we obtain Ga 2 m c <J^o? : /4 ~ 
14 eV, where m c is the carbon atom mass [27||. The condi- 
tion G' > then implies that Da < y/ Gahvppy w 9.8 eV. 
The experiment [3] suggests that valley splitting is larger 
than Zeeman splitting, which occurs in our model if 
Da > 6.3 eV, as we show below. It is quite difficult 
to accurately predict the constant D and we will there- 
fore only provide an order of magnitude estimate. The 
mechanism that we think gives the largest contribution 
results from the interaction of a single carbon atom with 
the SiC>2 substrate treated as a dielectric continuum [28j |. 
The non-retarded Lennard- Jones interaction energy of 
an atom at a distance r of a dielectric wall is given by 
S L j(r) pb -(er - l)(d 2 )/(e r + l)487re r 3 , where (d 2 ) is 
the atomic ground state expectation value of the squared 
electric dipole moment [18J|. The on-site energy change 
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resulting from the lattice distortion may be estimated as 

€ r - 1 (d 2 \ 



± M w E LJ (d ± r?) - S LJ (d ) « ± 



e r + 1 167reorfo 



t, (8) 



where the ± sign refers to sublattice A (+1) or B (—1) 
[jjj ]. do is the average distance between the graphene 
sheet and the substrate and we assumed that r\ <C do. 
For a carbon atom \J (d 2 ) ~ 4ea , where a is the Bohr 
radius, which gives Da ~ a(e r — l)e 2 ag/ (e r + l)7re dQ ~ 1 
to 14 eV depending on do ~ 0.1 to 0.2 nm. Therefore, 
the order of magnitude of the deformation potential Da 
is 5 eV. From now on, in order to match the experiment 
3, we take the plausible value Da — 7.8 eV, which gives 
G'a 2 « 4.2 eV. 

Minimizing E tot = E n=0 + E B<0 + Elastic as a function 
of the distortion 77, we obtain an on-site energy 



M = Drj = 
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and a condensation energy E tot — — (2— \v\)N,pM /2. The 
distortion is indeed very small, of order rj/a ~ 2.10 -5 x 
Sj_[T] when v w 0. This gives an n = valley splitting 
A v = 2M w 4.2K x (1 - |i/|/2)5j_[T], which for v w 
is larger than the Zeeman splitting = 3*Ms^tot ~ 
1.5KxB tot [T] The on site energy M is indeed much 
smaller than the cyclotron energy and can therefore be 
safely neglected in each n 7^ LL: M/hw c ~ 5.10~ 3 x 
yj B± [T] when 1/ ss 0. This means that the LL spectrum 
for n 7^ is approximately given by Eq. (JJ) - as in the 
case of no lattice distortion - and therefore e n ^ a oc y/B± 
in agreement with recent spectroscopic observations [ijj. 

Considering LL broadening due to disorder, the pre- 
ceding calculation for lattice distortion is modified at 
weak magnetic field, when the valley splitting A v is 
smaller than the LL width A; mp . For example, for rect- 
angular LL - the density of states being 4iV^ / A; mp inside 
a LL and zero otherwise - Eq. ((4|) is changed into 
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while Eq. (0 remains unchanged since it concerns to- 
tally filled LLs. The electronic energy gain is now pro- 
portional to rj 2 . Comparing this energy to the renor- 
malized elastic energy loss of Eq. (jTj) , we see that a dis- 
tortion only occurs if 2N < pD 2 / A; mp > N P G' , which al- 
ways happen at large enough magnetic field. This con- 
dition is precisely equivalent to requiring that the val- 
ley splitting A„ = 2M - given by Eq. © with v w 
- be larger than the LL width A; mp . This is satisfied 
if B ± > hG'A imp /3V3ea 2 D 2 ~ 7 T, where we used 
A; mp = 2r with r ~ 15 K the measured LL half-width 
at half-maximum Therefore, as soon as B± is larger 
than this threshold value, the lattice is distorted and the 
valley gap is larger than the LL width, which means that 
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Figure 1: Energy e of the first LL versus magnetic field B. 
The degeneracy in units of the flux number N$ appears on 
the levels. The cyclotron hu c , valley A„ and Zeeman A z gaps 
are also specified. At large B, the levels are tagged by the LL 
n, spin a and valley a indices: (n, a, a). 



one can use the results obtained in the preceding para- 
graph in the case of infinitely narrow LL. 

We now discuss the expected plateaus in the Hall con- 
ductivity a xy = ve 2 jh as a function of the filling fac- 
tor v oc Vg and the magnetic field. We consider a sys- 
tem at low temperature T <C Aj mp , where thermal ef- 
fects can be neglected, and assume broadened LL with 
a width A; mp ~ 30 K that we compare to the calcu- 
lated gaps: for typical magnetic fields, the cyclotron gap 
huj c w 420K x \/B~x [T] is the largest, then the valley 
gap is A v » 4.2K x (1 - |z/|/2)S x [T] and finally the 
Zeeman gap Az ~ 1.5K x £? tot [T] is the smallest, see 
Figure [TJ When the magnetic field is such that the cy- 
clotron gap becomes larger than Ai mp which occurs at 
B± r~j 5.10~ 3 T, one expects plateaus at v = ±(4|n| + 2). 
Then, when the valley gap (corrected by Zeeman split- 
ting) A v — Az becomes of order A; mp which occurs at 
B± r~j 11 T for i/RiO - thanks to our choice of Da - one 
expects a v — plateau. Finally, when the Zeeman gap 
reaches A; mp which occurs at B tot ~ 20 T, one expects 
plateaus at v — ±1 and v — An (n ^ 0). Because val- 
ley degeneracy is not lifted by the lattice distortion when 
n 7^ 0, plateaus are not expected at v — ±(2|n| + 1) 
(n 7^ 0). Experimentally, plateaus at v = ±2; ±6; ±10 
are observed at magnetic field ~ 9 T [l|, 0| and are at- 
tributed to the cyclotron gap, the v = plateau appears 
at 11 T, v = ±1 and ±4 are observed at B± > 17 T 
and the v = ±3; ±5 plateaus are not observed Q. This 
agrees qualitatively with our model and allows one to at- 
tribute the v = plateau to the n = valley gap and the 
v = ±1; ±4 plateaus to the Zeeman gap. 

The v = plateau, which occurs around zero gate volt- 
age, is worth considering from an edge states perspective 
[20(1 . We assume smooth edges on the sides of a sample of 
width W and take infinite mass confinement as bound- 
ary condition, following Ref. [2l|. The on-site energy 
(the "mass") is now position dependent in the y direction 
perpendicular to the edges: in the bulk, M{y) is con- 
stant and given by Eq. ((9]); on the edges y w ±W/2, it 
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smoothly rises to infinity in order to confine the electrons. 
Eq. |(2| and ([3} show that electronic states with positive 
(resp. negative) energy bend upward (resp. downward) 
in energy on the edges as M(y) — > oo. As > A 2 , 
the sign of the energy is given by that of the LL index n 
except for n — where it is given by the valley index a. 
Therefore when the Fermi level lie in the valley gap, there 
are no edge states, and the Hall conductivity a xy = 0, 
as expected. The absence of edge states (when v « 0) 
is a consequence of the valley splitting being larger than 
the Zeeman splitting, see Ref . @, 0| . At the same time 
the longitudinal conductivity a xx should be exponentially 
small (activated) because of the absence of current car- 
rying states both in the bulk and on the edges: therefore, 
one does not expect a wide zero in the longitudinal re- 
sistivity l/a xx , as for usual QH states, but rather in the 
conductivity a xx . Actually, the system should conduct 
as a very bad metal, which according to Mott's criterion 
implies a xx ~ e 2 //i, just as for graphene under zero mag- 
netic field [3, 0|. This point deserves further studies. In 
the experiment Q, when a v = plateau is observed in 
a xy at 25 T, the longitudinal resistance features a finite 
peak R xx ~ 40 kfi, corresponding to a resistivity of or- 
der l/a xx ~ 10 kfi of the same order as that measured 
at zero magnetic field l/a xx h/Ae 2 « 6.5 kfi [3,0]. 

In conclusion, we compare the predictions of our model 
to that of valley ferromagnetism @, SS First, we 

predict that valley degeneracy is not lifted in n ^ LL, 
whereas valley ferromagnetism lifts this degeneracy. This 
results in the absence of the v = ±(2|n| + 1) plateaus, 
with n/0. Second, the valley gap is proportional to 
the perpendicular magnetic field, whereas the n = 
skyrmion gap relevant for ferromagnetism A^y ~ e 2 /cIb 
scales as y/B± 0, HJ: this should be seen in activa- 
tion gaps measurements Q. In addition, using the co- 
incidence method with a tilted magnetic field [3], one 
should be able to distinguish the different gaps through 
their dependence in the perpendicular or total magnetic 
field. The gate voltage dependence of the valley gap 
A v oc (1— |f |/2) could be detected spectroscopically [1131 ] - 
Third, if lattice distortion indeed occurs it should be di- 
rectly seen. It might be detected using synchrotron X-ray 
diffraction at grazing incidence or scanning tunneling mi- 
croscopy at magnetic fields ~ 10 T and low temperature. 
Fourth, the lattice distortion and its consequences should 
vanish if the graphene sheet is placed in a symmetric di- 
electric environment. In the end, we provide what we 
think is a plausible mechanism for lifting valley degener- 
acy. Whether lattice distortion indeed occurs remains to 
be checked experimentally. 
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